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O ' ABSTRACT 

m \ 
O ■ 

\ The geometric Baade-Wesselink method is one of the most promising tech- 

niques for obtaining a better calibration of the Cepheid period-luminosity relation 
by means of interferometric measurements of accurate diameters. In this paper we 

\ present new wavelength- and phase-dependent limb darkening corrections based 

on our time- dependent hydrodynamic models of the classical Cepheid C Gem. 
. £h ! We show that a model simulation of a Cepheid atmosphere, taking into account 

the hydrodynamic effects associated with the pulsation, shows strong departures 
from the limb darkening otherwise predicted by a static model. For most of 
its pulsational cycle the hydrodynamic model predicts a larger limb darkening 
then the equivalent static model. The hydrodynamics affects the limb darkening 
mainly at UV and optical wavelengths. Most of these effects evolve slowly as 
the star pulsates, but there are phases, associated with shocks propagating into 
the photosphere, in which significant changes in the limb darkening take place 
on time-scales of the order of less than a day. We assess the implication of our 
model LD corrections fitting the geometric Baade-Wesselink distance of ( Gem 
for the available near-IR PTI data. We discuss the effects of our model limb 
darkening on the best fit result, and analyze the requirements needed to test the 
time-dependence of the limb darkening with future interferometric measurements. 



1 We are saddened to report that our colleague, Costas Papaliolios, died 2002 June 6. 



- 2 - 



Subject headings: Cepheids — stars: atmospheres — stars: oscillations — tech- 
niques: interferometric — stars: individual (£ Gem) 

1. Introduction 

With the recent advances in interferometry, the Baade-Wesselink (BW) method (Baade 
1926; Wesselink 1946) has become one of the most promising techniques to obtain inde- 
pendent distance measurements of pulsating stars. A "geometric" variation of this classical 
method (Sasselov & Karovska 1994) can in principle be used to derive the distances of nearby 
classical Cepheids from the variations, as they pulsate, of their angular diameter and radial 
velocity. The potential accuracy of this method can improve the calibration of the Period- 
Luminosity relation of Classical Cepheids (Leavitt 1906), which is a fundamental step in the 
cosmological distance ladder. 

The geometric BW method has been recently applied on two Cepheids in the northern 
hemisphere, ( Gem (Lane et al. 2000) and rj Aql (Lane et al. 2002). Both stars were 
observed at near-IR wavelength (H-band) with the Palomar Testbed Interferometer (PTI). 
The uniform brightness disk diameter was measured at several epochs for both stars, with 
enough precision to detect the changes in the stellar radius due to the pulsation. The 
measured angular diameters were then fitted with an appropriate model derived from the 
radial velocities of the two stars, obtaining the distance with an estimated accuracy of ~ 10%. 

As explained in detail by Lane et al. (2002), the main source of uncertainty in the 
measurements is related to (1) the conversion of photospheric line velocities into radial motion 
(the so-called projection factor, or p-factor), and (2) the estimate of the limb darkening 
(LD). These two quantities can be derived by modeling the Cepheid atmosphere, taking 
into account the spherical geometry of the star and the time-dependent hydrodynamics of 
the pulsations. The correct way to determine the p-factor for pulsating Cepheids have been 
described in detail in Sabbey et al. (1995) and Krockenberger et al. (1997). More recently, 
we have presented a new method for computing accurate time- and wavelength-dependent 
center-to-limb brightness distributions for classical Cepheids (Marengo et al. 2002, hereafter 
paper I). 

The model described in paper I provides a significant improvement of the limb darkening 
coefficients with respect to the tabulated values currently used to analyze interferometric 
data. These tables (Parsons 1971; Manduca 1979; Kurucz 1993a; Claret et al. 1995) are 
based on hydrostatic model atmospheres of non-pulsating yellow supergiants, having similar 
T cff and log g as classical Cepheids. Our models have the advantage of being specific for each 
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simulated Cepheid, and provide the appropriate limb darkening for any pulsational phase at 
arbitrary wavelength. As shown in paper I, the limb darkening predicted by our models is 
significantly different from the one expected for a static yellow supergiant. 

In this paper we analyze how the inclusion of hydrodynamic effects in our Cepheid 
models can affect the interferometric distance determination with the geometric BW method. 
We first compute the limb darkening corrections as a function of wavelength and pulsational 
phase (section 2), starting from the center-to-limb intensity profiles presented in paper I 
for the classical Cepheid ( Gem. Then, we derive a new estimate for the geometric BW 
distance of this star, correcting the PTI data presented by Lane et al. (2002) with our phase 
and wavelength dependent LD corrections (section 3). The accuracy of the best fit values, 
and the importance of various error sources, are described in detail. We then conclude by 
analyzing the extent of the corrections induced by our models in the geometric BW distances, 
and the level of accuracy that will be required by present and future interferometers to detect 
the hydrodynamic effects predicted by our models. 



2. Limb Darkening Corrections 

The models described in paper I consist of a series of synthetic atmospheres computed 
with second-order one-dimensional hydrodynamic calculations, performed assuming a spher- 
ical symmetry. Each model simulates the atmosphere of a Cepheid at a certain pulsational 
phase. In paper I we have described the general procedure to obtain such models, and have 
shown the specific results for the classical Cepheid ( Gem. 

The pulsational period of ( Gem (of approximately 10.15 days) was covered by a total 
of 49 individual models (the time-step in our model grid is determined by the convergence 
criteria in the hydrodynamic simulation). The time resolution of our model sequence is thus 
~ 0.2 days. For each model in the sequence, we have computed realistic spectral intensity 
distributions. This step was done by approximating the dynamic models with a plane-parallel 
atmosphere. The radiative transfer problem for the static atmosphere was then solved in 
LTE conditions using the ATLAS code (Kurucz 1970, 1979, 1993b) and its opacity library. 
The end result of this procedure is a set of 49 spectral intensity distributions describing the 
( Gem spectrum as it changes while the star pulsates. 

In paper I we showed how this procedure allows us to compute accurate limb darkening 
profiles for ( Gem, which are (pulsational) phase and wavelength dependent. We proceed 
here to the next step, which is to describe how these models can be used in interferometry. 

Interferometric measurements are usually expressed in terms of uniform intensity disk 
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(UD) or limb darkened diameters. The UD diameters are derived by fitting the normalized 
fringe visibilities (squared) V 2 with a uniform disk model, which assumes that the star is a 
disk of uniform brightness: 



where 6ud is the UD angular diameter and p the spatial frequency of the measurement. Since 
a uniform disk model is derived assuming that the source is a disk of uniform brightness, it 
does not depend on the wavelength. Real stars, however, are limb darkened and have limb- 
to-center brightness distributions which are a function of wavelength, as shown in paper I. 

For any practical purpose, we can consider the physical radius of the stellar photosphere 
to be a well defined quantity. The contribution functions of observable spectral features in 
the stellar spectra does peak at different heights in the atmosphere, but the difference is very 
small compared to the full diameter of the star. For this reason, we can assume that the 
angular diameter of a Cepheid is the same at all wavelengths, after being corrected for LD. 
This diameter, which we call "limb darkened diameter", since it takes into account 

the fact that the star has LD, will thus only depend on the pulsational phase of the star, 
and not on A. Fitting fringe visibilities with a UD model, on the other hand, will produce 
a different 9 UD according to the wavelength, because of the inability of UD models to take 
into account the spectral properties of the LD. For this reason the UD angular diameter 
measured from fringe visibilities of a LD star will instead depend on both the pulsational 
phase and the wavelength: 9ud(\ </>)■ 

To obtain the real angular diameter of a star from the fringe visibilities, one can either fit 
the V 2 data with an appropriate limb darkening model, or correct the UD measurement with 
a specific wavelength dependent limb darkening correction. The first approach is certainly 
preferable, because preserves the wavelength dependence of the measured visibilities, and 
thus allow a test the limb darkening models. There are cases, however, in which the original 
visibilities are not available, as for example when mixing heterogeneous data from different 
interferometers. 

In this section we present accurate LD corrections derived using our models. They are 
based on hydrodynamic models and therefore they are phase and wavelength dependent. 
Even though the models are specific for £ Gem, we use them here as an example to de- 
scribe the more general case of LD corrections for the generic classical Cepheids, and their 
implications for interferometry and the BW method. 




(1) 



The conversion factor between UD and limb darkened diameters is usually defined as 
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k(\,<f>) = 9ud(^, <t>) / Q ld{4>) ■ This factor is in many cases approximated by the model of a 
star having similar spectral type as the source of interest (in the case of Cepheids, a yellow 
supergiant, as described in Claret et al. 1995). Such approximations do not take into account 
the phase dependence for variable sources, and may be incorrect in their dependence on A, 
since the atmosphere of a Cepheid cannot be properly approximated with a static supergiant 
(see discussion in paper I). 

The k(\, 0) corrections discussed in the following sections are derived with the following 
procedure. At each phase in the model grid, and for a set of wavelengths in the optical and 
near-IR range, the center-to-limb profiles computed in paper I are converted into fringe 
visibilities by applying the Hankel transform (see e.g. Koechlin 1988): 



where p is the spatial frequency of the interferometric visibility, w = (6ld/2) • sin(a) is 
the projection of the stellar angular radius on the disk and I u (w, A, 0) is the model stellar 
intensity spectrum at the projected radius, computed for the given pulsational phase 0. 

The simulated Vl D visibilities are then fitted with an UD V^ D model, in order to derive 
the UD diameter Oud which the interferometer would have measured, at each wavelength, 
for the model source. The limb darkening correction A;(A, 0) is the ratio between the best 
fit 9ur> and the value of the angular diameter Q^d = 1 mas that we set to give a physical 
dimension to the simulated visibilities in the Hankel transform. The simulated visibilities 
were computed over a spatial frequency range up to the first minimum, to match the typical 
conditions encountered when fitting real data with currently available interferometers. We 
have however tried several other combinations of model 9 LD and spatial frequency ranges 
(to simulate different interferometer baselines), confirming that the resulting A;(A,0) are not 
very sensitive to these parameters, as long as the simulated star is at least partially resolved. 



An example of the wavelength dependence of the LD correction A;(A) is shown in Fig- 
ure 1. The figure shows the LD correction derived from our hydrodynamic simulations (thick 
line). The thin line is the LD correction obtained from hydrostatic atmospheres (Kurucz 
models) having the T e ff(0) and logg(0) measured as a function of the pulsational phase by 
Krockenberger et al. (1997). Note that the T e g of the dynamic and static models, at each 
pulsational phase, is the same, since the observational effective temperature was used as an 





2.1. Spectral Properties of Limb Darkening 



- 6- 



input parameter to compute the dynamic model. The log g is instead different in the two 
cases, as a consequence of the procedure followed to solve the radiative transfer in the hy- 
drodynamic case, and due to the specific definition of the gravity terms in the hydrodynamic 
equations. 

The k correction is strictly related to the spectral properties of the source. The limb 
darkening is higher (lower k) at UV and visible wavelengths, and converges toward unity (no 
limb darkening) with increasing wavelengths toward the infrared. Spectral lines are less limb 
darkened than the continuum: having higher optical depth, they probe upper atmospheric 
layers, where the temperature gradient is lower. Among the spectral features that are visible 
in our low spectral resolution computations, is the Balmer jump in the UV, the Ca H&K 
doublet (A ~ 395 nm), H/3 (A ~ 486 nm) and Ha (A ~ 656 nm). 

The first panel in Figure 1 shows the limb darkening correction at minimum radius. At 
this phase, the Cepheid atmosphere is most compressed, giving rise to a steeper temperature 
gradient. This is responsible for an increase of the limb darkening at all wavelengths. The LD 
correction for the hydrostatic atmosphere, on the other hand, does not show this same effect, 
having a flatter limb darkening despite the fact that T cS is the same in the two models. The 
difference is as much as Ak ~ 0.01 at visible wavelengths, and less in the near-IR. Accurate 
measurements of the LD with optical interferometers should then be able to verify this effect. 
The difference between the static and dynamic models is even larger at the wavelengths of 
the main spectral features, which appear stronger in the more compressed hydrodynamic 
atmosphere. 

In the next section we will show that in most cases the effects of the hydrodynamics 
on the atmospheric structure result in a larger LD. There are however phases in which the 
free expansion of the atmosphere results in a quasi-static structure. In these phases (before 
and after maximum luminosity), the hydrodynamics is less important and the two k(X) are 
virtually indistinguishable. 

As shown in the bottom panel, however, the full force of a shockwave crossing the 
photosphere results in high excitation states and local expansion. Even though at this phase 
£ Gem is contracting (this happens one day before minimum radius), the energy deposited by 
the shock in the region where the visible photons are created generates a lower temperature 
gradient. As a consequence, the LD of the hydrodynamic simulation is lower than in the 
static case. Note, however, that this effect is mostly appreciable at visible wavelengths. 
An optical interferometer should detect at this phase a decrease of the limb darkening of 
Ak ~ 0.01. This phase is very brief and, at least in the case of ( Gem, the effects of the 
shock are already dissipated after less than 20 hours, when the LD is maximum again as the 
star approaches minimum radius. This timescale is related to the shock propagation speed, 
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which is very well constrained by observations (see Figure 8 of Sasselov & Lester 1994). 

2.2. Limb Darkening and Pulsational Phase 

By convolving the k(X, <fi) models with the filter passbands used by interferometers, we 
can study the variations of the LD correction as a function of the pulsational phase. For a 
convenient comparison with observations, the best choice is to adopt the phase as defined by 
the optical lightcurve, where the zero phase, 4>l = 0, coincide with maximum luminosity. As 
explained in paper I, however, our dynamic simulations are computed as a function of the 
phase <f)y based on the pulsational velocity, in which zero phase is at minimum radius. To 
convert our LD models to the lightcurve phases we have used the phase shift A0 = <Pl ~ <Pv 
between maximum luminosity and minimum radius computed by Bersier et al. (1994a,b). 

Figure 2 shows the limb darkening corrections plotted as a function of the visible 
lightcurve phase for two wavelengths. The top panel shows k(<f>) for the near-IR H filter 
passband used by the PTI interferometer. The bottom panel shows k(<t>) for the bluest avail- 
able channel (A ~ 570 nm) of the Navy Prototype Interferometer (NPOI) in the optical. 
The curves have been filtered to remove numerical noise with an adaptive gaussian kernel, 
and the 1 rms error bands in our LD corrections associated to the numerical uncertainty are 
indicated by two thin lines bracketing each LD curve. 

The plots show that k(<f>) is roughly constant for most phases, with a slow increase as the 
star expands, as a consequence of the decreasing temperature gradient in the photosphere. 
Coincident with the passage of a shockwave through the photosphere (at <p L ~ 0.6), k{4>) 
shows a sudden rise in the optical, indicative of a sharp decrease in the limb darkening. 
The energy deposited by the shock is responsible for this effect. It increases the excitation 
in the photospheric layers, and thus flattens the temperature gradient. This extra energy 
dissipates in a short timescale, after which the atmosphere resumes its normal state. At 
(f) L ~ 0.7 ( Gem is close to minimum radius, where the temperature gradient is higher (due 
to the compression of the atmosphere) and thus the star appears more limb darkened. 

The sharp increase in k(<p) observed in the optical, at the time of the shockwave, is 
not predicted for the H band. This is because in this wavelength range the effects of the 
shock are less pronounced. This might be due to the fact that the emergent spectrum in the 
H-band originates deepest in the atmosphere as a result of the broad minimum of the H 
opacity at 1.6 /iin, and the shocks increasingly steepen and disturb the atmosphere as they 
propagate down the density gradient (up into the atmosphere). 

It is important to understand how much of this behavior in the LD curves is due to 
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hydrodynamic effects, or just to the changes in T c g and logg during the stellar pulsation. 
Solid lines in Figure 3 show our hydrodynamic simulations for two wavelengths. The dotted 
lines show the LD correction derived for hydrostatic atmospheres having at each phase the 
T e gf and log g measured by Krockenberger et al. (1997), which we used as starting points in our 
models. The overall value of k{<p) is similar, at each wavelength, between the hydrodynamic 
and hydrostatic atmospheres. The phase dependence is instead completely different. Figure 4 
reveals that static model LD corrections closely follow the change of T eS with phase. In a 
hydrostatic atmosphere a higher T cS is responsible for a flatter photospheric temperature 
gradient, and thus for less limb darkening (i.e. higher k). In a hydrodynamic atmosphere the 
temperature gradient is determined predominantly by the pulsation dynamics, and therefore 
the phase dependency of the LD is different. In a hydrostatic atmosphere, gravity seems 
to have little or no effect determining the changes in the limb darkening. When computing 
the spectral energy distribution of our dynamic atmosphere, we constrain the model T c g 
to the observed values, leaving \ogg as a free parameter to be fitted with a grid of static 
atmospheres (see paper I). The rationale for this choice was that T e g has a dominant role in 
determining the spectral properties of the atmosphere. This result confirms the validity of 
our choice. 

To test the assumption that logg does plays a minor role, we have also computed the 
LD correction for the static models having both T c g and log g identical to the hydrodynamic 
simulations. The results are very similar. Despite the significantly different logg with 
respect to the observed values, k{4>) still closely follows the variations of T cS . Only just after 
minimum radius (when T eff is maximum and the best fit logg is minimum), logg plays some 
role in lowering the LD correction. This confirms that, except when logg is very small (~ 1, 
in this case), the effective temperature is the decisive parameter for the LD in a hydrostatic 
atmosphere. 

In a hydrodynamic atmosphere, however, T eS is not the dominant parameter. As shown 
in Figure 3, when hydrodynamic effects are taken into account, the relation between T e g 
and the LD breakes. Shocks are the source of the most dramatic effects in k(<f>), but even 
when they are absent, the LD corrections are strikingly different from the ones predicted 
by hydrostatic atmospheres. This is because they are not determined by T eS , but by the 
time-dependent structure of the expanding/contracting atmosphere. Given that our hydro- 
dynamic models generally have a steeper temperature gradient than a similar hydrostatic 
atmosphere, the resulting LD is larger for most phases. 

Figure 3 also shows the LD corrections used in Lane et al. (2002) for ( Gem, which has 
the constant value of k — 0.96 ±0.01. This correction has been derived from tables published 
by Claret et al. (1995), based on hydrostatic models. Note that the average value of our LD 
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correction is outside the error bars of the tabulated value. The reason of this discrepancy 
is that our model closely follows the T eS and log g derived for ( Gem from spectroscopic 
observations, while the Claret et al. (1995) values are computed for a generic grid of yellow 
supegiants matching the Cepheids spectral type. 

Note finally that, as expected, the LD correction is much smaller at near-IR wavelengths, 
and larger toward the blue. At bluer wavelengths the amplitude of the k variations with 
phase (especially at the time of the shockwave) is also much larger. Therefore, to test 
observationally the details of the hydrodynamic effects on limb darkening, interferometers 
operating at short wavelengths are favored. 



The LD deviations introduced by the hydrodynamics are relatively small (1% in the 
optical, and less in the infrared), and are thus usually ignored when treating interferometric 
measurements. In order to obtain accurate distances with the BW geometric method, how- 
ever, higher levels of accuracy are needed. Under such stringent requirements, the time and 
wavelength dependent hydrodynamic effects play an important role, and should be taken 
into account. 

A detailed description of the problems involved with the application of the geometric 
BW method to pulsating Cepheids is given in Sasselov & Karovska (1994). This method 
allows the determination of the distance and the average radius of a pulsating star from both 
its angular diameter and radial velocity at several phases. This is done by means of a x 2 fit 
of the following function: 



where 6 = 2R Q /D is the average angular radius of the star, and ARi its radial displacement 
at the time of the interferometric measurement 0j. The variation of the stellar radius AR((f) L ) 
is derived, as a function of the lightcurve phase 0l, by integrating the pulsational velocity 
over time: 



3. Effects of the LD corrections on the geometric BW method 




(3) 




(4) 
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where v r (<f>) is the radial velocity, which is corrected by the systemic velocity 7 and the 
p-factor to yield the true pulsational velocity The systemic velocity can be derived by 
requiring the conservation of the radius over one period (see paper I). Appropriate p-factors 
have been computed for pulsating Cepheids such as ( Gem; we use here the value of 1.43 
(the same adopted by Lane et al. 2002), derived from the hydrodynamic model we used to 
compute our LD profiles, which was published in Sabbey et al. (1995). 

The fit in equation 3 can be solved analytically. The best fit Rq and D can be derived 
by minimizing the x 2 relation solving for the two unknown parameters. Note that we do not 
make use of the color curve of the Cepheid, which can be applied to further constrain the 
geometric BW solution. A detailed description of this procedure will be given in a separate 
paper. 

The LD corrections k((f>) computed in the previous sections enter the fit by allowing the 

(i) 

conversion of the UD diameters 6 UD obtained by the interferometer into true LD diameters 
of the star: 0; = 

measurement, computed for the wavelength of the observation. 



0u' D /k((f>i), where k((f>i) is the LD correction at the phase of the i-th 



The terms Oi in the fit \ 2 relation are the errors associated with each data point. They 
are the geometric sum of the individual error sources (in terms of angular diameters) for 
each observation: 



a. 



+ 



o. 



(AR/D)' 



(7. 



(fc) 



(5) 



where a\ is the error of the interferometric measurements, a\ the error related to the 
radial displacement and af^ the error in the LD correction. According to Lane et al. (2002), 
the errors in the PTI H-band data vary between 0.01 and 0.06 mas. The errors on the radial 
displacements are due to the uncertainty in the p-factor (amounting to ~ 4%, according to 
Sabbey et al. 1995), and to the estimated measurement errors in the radial velocity, adding 
a further 2% according to Bersier et al. (1994b). The errors in our LD corrections, barring 
systematic errors and based only on the numeric uncertainty in our model, are of the order 
of ±0.02% (see Figure 2), which is a negligible contribution with respect to the other error 
sources, and is thus ignored. 

The best fit results are shown in Table 1 (col. [2]) and Figure 5. Column (3) in Table 1 
shows, for comparison, the best fit parameters obtained with the same fixed LD correction 
used by Lane et al. (2002) of k ~ 0.96. The difference in the two best fit values of the 
geometric BW distance is 11 pc, which is less than ~ 3%. Even though this difference is 
less than 1/3 of the error bars, it is significant: Figure 5 shows that the error regions for our 
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best fit, the UD fit and the fixed k fit are mutually exclusive. This is again a consequence 
of the high level of precision that the available interferometric data already allow in the 
determination of the average angular radius 6 . 

Note that the error regions shown in Figure 5 do not take into account the uncertainty 
in the LD correction, which we ignored. The uncertainty quoted for the tabulated LD 
correction used by Lane et al. (2002) is much larger than our numeric error Ak ~ 0.02, and 
would have resulted in a much larger error region, including all the narrow ovals in Figure 5. 
The main point of this paper, however, is to discuss the consequences of using model derived 
LD corrections, and Figure 5 shows that in the case of ( Gem using generic tabulated LD 
corrections can lead to an error of a few percent with respect to the best fit distance. This 
discrepancy is larger than the error regions determined by all the other error sources, and 
can thus in principle be tested observationally by directly measuring the LD with longer 
baseline interferometers. 

The difference in the results of the geometric BW fit are entirely due to the LD cor- 
rection, more precisely to its average value k computed over the pulsational cycle. In the 
near-IR, and especially in the H band, the phase dependency of k(\, 0) is relatively small 
(less than 2%, compared with the error bars > 5 — 10% in the measured angular diameters). 
This means that, with the current error bars in the interferometric data, we are still not sen- 
sitive to the variations in the LD corrections induced by the hydrodynamics. An accuracy 
of the order of 0.2% in the 6i is required to be sensitive to such effects in the H band. Note 
that a more favorable situation is met in the visible, where the temporal variations of k((j)) 
are larger. Given the limitation of the available data, however, the main contribution of our 
hydrodynamic simulations, at least for now, is not in the detailed dependence of /c(0), but in 
setting the right level of the LD correction for the wavelength of the observation. Contrary to 
the tabulated values of the LD correction, which are computed for generic yellow supergiants 
having the same spectral type of Cepheid stars, our corrections are specific for the modeled 
stars, as they follow the T e s and log g derived from spectroscopic observations for each star. 
This is not a trivial matter, as the discrepancy with the tabulated value (AA; ~ 0.02 in the 
H band) is the largest error source in the determination of the geometric BW distance. 

Finally, we discuss the issue of the validity of the best fit x 2 as guide to which LD better 
represent the data. Is the best fit x 2 an indicator of the best value for the LD correction? If 
this is the case, then we should conclude that ( Gem is not limb darkened, since the best fit 
X 2 computed with this same procedure with the UD diameters would be x 2 — 27.5. However, 
this is not the case. 

To explain this, one should first remember that the time dependent hydrodynamic effects 
are too small to be observed with the available PTI data. While doing the BW fit, we could 
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thus consider the LD correction as a constant value k equal to its average over the pulsational 
period. In the three fits shown in Figure 5, the value of k would be equal to 1 for the UD 
data, ~ 0.979 for our hydrodynamic calculations and 0.96 for Claret et al. 1995 (used by 
Lane et al. 2002). This coefficient enters the x 2 equation by dividing the UD 9jj D : 



* 2 = £ 

i 

A smaller value of k (larger LD), will result in a larger best fit mean angular diameter 
O , which will thus scale as 1/k. The net effect of this scaling is that the best fit x 2 will 
scale as 1/k. This means that a fit made with larger LD (as long as the phase dependency 
of k is irrelevant) will have a larger best fit x 2 ■ The UD fit will thus appear always to be the 
best. 

The conclusion is that, until the accuracy of the measured angular diameters becomes 
good enough to appreciate the changes in LD due to the pulsation, the geometric BW fit 
is not a good tool to test the LD models. On the other hand, a good model for the LD is 
absolutely necessary to obtain a reliable value of the geometric BW distance, since the best 
fit angular diameter scales linearly as 1/k. 

Note, finally, that a different fitting strategy may allow a test our LD models even with 
present day interferometric data. Our insensitivity to the LD is due to the fact that we are 
fitting derived data (the angular diameters 0j) which have been previously obtained from the 
fringe visibilities assuming a uniform disk model. This initial fit has destroyed the sensitivity 
of the original data to the stellar LD. The correct way to test the LD is to perform x 2 fits on 
the original visibilities. Any test for the absolute value of LD predicted by our models should 
thus be made using the visibility data taken at different projected baselines and wavelengths. 
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4. Conclusions 

The geometric BW method is a powerful tool to derive the distances of pulsating stars. 
The detailed analysis of the method, and our discussion of its main current uncertainties 
described in the previous sections, show that a special care should be taken when assessing 
the accuracy of the results. 

Despite the increasing quality of the interferometric data, the error bars in the individual 
measurements are the largest contributors to the final errors in the geometric BW distances 
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and average radii. After the instrumental errors, however, follows the uncertainty in the 
limb darkening correction. The importance of LD will grow in the near future, given the fast 
pace at which the available interferometers are improving their accuracy, and with the new 
long baseline, large aperture interferometers which are becoming operational. To address 
this issue, we have presented in this paper a procedure to compute accurate LD corrections 
for pulsating Cepheids, which are based on time- and wavelength-dependent hydrodynamic 
models. 

We show that our k(X, <fi) strongly differ from the equivalent corrections computed from 
hydrostatic atmospheres, even in absence of "strong" hydrodynamic effects like shocks. The 
main consequence of the hydrodynamic terms in our stellar atmospheric models is an increase 
in the limb darkening, because of a generally higher temperature gradient in the photosphere. 
The situation briefly reverses in presence of Shockwaves, due to the energy deposited by the 
shocks. 

The magnitude of these effects depends on the wavelength. The effect is more significant 
in the visible spectrum (Ak ~ 0.015) and it is smaller in the infrared (Ak ~ 0.002). Cur- 
rent interferometers still cannot test these time-dependent variations in the limb darkening, 
which will however become important in the near future, especially when interferometric 
observations in spectral lines will become feasible. The phases of rapid variation in the LD 
associated to the propagation of the shockwave in the photosphere, when significant changes 
in the LD occur in timescales of hours, appears particularly appealing for an observational 
test of our models. 

Even though the time-dependent variations in the LD are not currently measurable 
because of the required visible accuracy, our models can already provide an average value 
of the wavelength-dependent LD correction k(X) which is significantly different from the 
tabulated values currently used. In the case of ( Gem, our limb darkening corrections induce 
a ~ 3% change in the best fit value of the BW distance derived from H-band PTI data. 
Even more important is that, as shown in section 3, the la error regions around the best 
fit distances derived with our and other LD corrections are mutually exclusive, opening the 
possibility of an independent test of our models when longer baselines reaching the first 
minimum in the visibilities will allow a direct test. A preliminary test of our models will also 
be possible by directly fitting the observed visibilities with our model visibilities, instead of 
using LD corrections of UD best fit diameters. 

As the data for other Cepheids will become available, it will be important to have a large 
library of models specifically computed for each source. As the models are very dependent on 
the pulsational characteristics of each star, upon which the hydrodynamic model is build, a 
"generic" parametric limb darkening correction for all stars is not possible, as any individual 
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model cannot be extended to be used for a Cepheid with a different pulsational engine. This 
is the reason why tables of LD corrections as the ones produced by Claret et al. (1995) 
cannot reproduce the detailed changes in the LD which are required to apply the geometric 
BW method to classical Cepheids. 

We wish to thank the anonymous referee and the editor for the comments and sugges- 
tions that helped us to improve this paper. This work was partially supported by NFS grant 
AST 98-76734. M.K. is a member of the Chandra Science Center, which is operated under 
contract NAS8-39073, and is partially supported by NASA. 
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TABLE 1 

BW BEST FIT OF ( GEM PTI H-BAND LD DATA 

Fit Parameter k(<f>) (model) k = 0.96 

0) (2) (3) 

R [R ] 67.0l|;J 66.2±|;| 

L> [pc] 372lg 36lt£ 

O [mas] 1.665±0.007 1.698±0.007 

X 2 28.1 29.5 
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Fig. 1. — Wavelength dependence of the LD correction k(X) for our ( Gem model. From top 
to bottom, the models are shown at minimum radius, at a quasi-static phase after maximum 
luminosity, and at the time in which a shockwave is crossing the atmosphere. Solid thick line 
is k(X) from our hydrodynamic model, while the thin lines are derived from a hydrostatic 
model having the T e g and logg determined from observations (Krockenberger et al. 1997). 

Fig. 2. — Pulsational phase dependence of the LD correction k(<f) for our ( Gem model. The 
LD correction is shown for the near-IR H band (top), and at visible wavelength (570 nm, 
bottom). The thin lines are the 1 rms numerical uncertainties in our simulations. 

Fig. 3. — Phase dependent fc(0z,) for ( Gem in the PTI H band, and in the bluest 570 nm 
channel of NPOI. Solid lines are the LD corrections computed with our hydrodynamic model; 
dotted lines are the equivalent corrections for an hydrostatic atmosphere having T e g and logg 
from Krockenberger et al. (1997). The dashed line is the value computed by Claret et al. 
(1995) and used in Lane et al. (2002). 

Fig. 4. — Effective temperature, as a function of lightcurve phase, used to compute the 
hydrodynamic model of ( Gem. The dependence of T e g from the pulsational phase is derived 
from measurements by Krockenberger et al. (1997). 

Fig. 5. — Best fit parameters and error regions for ( Gem PTI H-band data. Top curve is 
the UD data, middle curve is our best fit for model LD data and the bottom curve is the 
result obtained using a fixed LD correction of k ~ 0.96 as in Lane et al. (2002). The inner 
error region is the 68% confidence level of the fit (Icr), while the outer is the 90% confidence 
level. 
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